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MULTIVARIABLE WAVELETS 



EVA CURRY 

Abstract. We investigate the connection between radix representations for 
Z n and self-affine tilings of R n . We apply our results to show that Haar-like 
multivariable wavelets exist for all dilation matrices that are sufBcie ntly large. 



1. Introduction 

We investigate the connections between radix representations for Z™, self-affine 
tilings of K™ , and Haar-like scaling functions for multiresolution analyses and asso- 
ciated wavelet sets. 

In a separate paper, we investigate the idea, also introduced by Jeong [5], of 
radix representations for vectors in Z", or general point lattices T — M(Z n ) (M a 
nondegenerate n x n matrix) 1 . We wish to consider expanding matrices which 
preserve T. Without loss of generality, we may assume r — Z™. A matrix which 
preserves Z n must have integer entries. 

Definition 1. A dilation matrix for Z™ is annxn matrix A with integer entries, 
all of whose eigenvalues A satisfy |A| > 1. 

Note that for a dilation matrix A, q :— \ detA\ is an integer, with q > 1. Then 
l/ 1 /A{X n ) has nontrivial cokernel. Let D be a complete set of coset representatives 
of Z n / J 4(Z n ). We call the elemen ts of D digits. 

We may associate a sequence of digits with each x E Z™ by the Euclidean algo- 
rithm, as follows. Each x € Z n is in a unique coset of Z n /A(2, n ), thus there exist 
unique x\ G Z" and ro G D such that 

x = Axi + tq. 

Similarly, for each Xj, j > 1, there exist unique Xj+i € Z™ and rj G I? such that 

Xj = Axj + i + r,j. 

Formally, we write 

oo 

x - ^2 AJr i- 

If there exists a nonnegative integer N such that rj — for all j > 2V, then the 
Euclidean algorithm terminates and we say that x has a radix representation with 
radix A and digit set D. 
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Definition 2. Let A be a dilation matrix. We say that the matrix A yields a radix 
representation with digit set D if for every ieZ" there exists a nonnegative integer 
./V = N(x) and a sequence of digits do, di, ■ . ■ , d^ in D such that 

N 
3=0 

That is, a dilation matrix A yields a radix representation with digit set D if 
every i£Z™ has a radix representation with radix A and digit set D. 
Let A be a dilation matrix, and define 

fi = min {cr : er a singular value of ^4}. 

Let F be a fundamental domain for Z™, centered at the origin, 

i r 

2' 2, 

In pp, we give the following two results about radix representations. 

Theorem 3. Let A be an n x n dilation matrix. If fj, > 2 £/ien A yields a radix 
representation of7L n with digit set D — A(F) H Z ra . 

Corollary 4. .For every dilation matrix A, there exists a positive integer /3 > 1 
smc/i £/ia£ A* 9 yields a radix representation with digit set Dp — A^(F) n Z". 

2. Radix Representations and Tilings 

Radix representations are closely related to sclf-afhnc tilings of R n . 

Definition 5. A measurable set Q C M. n gives a self-affine tiling of R™ under 
translation by Z™ if 

(1) Ufeez- {Q + A;) = M n , and the intersection (Q + fci) fl (Q + fc 2 ) has measure 
zero for any two distinct k\ , ki £ Z n (tiling) ; and 

(2) there is a collection of q = | det j4| vectors fci, . . . , k q € Z" that are distinct 
coset representatives of Z n /yl(Z n ) such that 

~ U? =1 (Q + fci) (self-affine). 

Set 

OO 

T = T(A, Z?) := {£ G M™ : £ = ^ A^'^} 

j=i 

with the digits G F> for some digit set D. One can easily check that T is a self- 
affine set. We would like to be able to think of the elements of T as the fractional 
parts of vectors in W 1 in the same way that the fractional par ts of real numbers lie 
in [0, 1]. This is an accurate interpretation if T is congruent to R"/Z™. Thought of 
another way, we would like T to tile M™ under translation by Z™. 

Theorem 6. Let A be a dilation matrix, and let D be a digit set for A. Then A 
yields a radix representation with digit set D if and only if the set T(A, D) tiles W l 
under translation by Z™ and the origin is in the interior of T . 

We split the proof of this theorem into a few lemmas. 

Lemma 7. Let A be a dilation matrix, and let D be a digit set for A. If A yields 
a radix representation with a digit set D, then the set T(A, D) tiles K™ under 
translation by Z n . 
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Proof. Applying Proposition 5.19 from [12] the dilation matrix A, digit set D, and 
corresponding set T = T(A, D) satisfy 

(1) T is a compact subset of R"; 

(2) A(T) = U deD (T + d); 

(3) U xeZn (T + x) = R"; and 

(4) T contains an open set. 

To show that T tiles M™ under translation by Z™, we must show that 

m((T + x) f](T + y)) = for all x ^ y, x, y G Z" 

(where m(-) denotes Lebesgue measure). We extend an idea of Lagarias and Wang 
([6], p.31) to show that m((T + x) n (T + y)) = for any distinct G Z" for 
which there is a radix representation with radix A and digit set D. 

Note that for any subset Q of R n , m(A(<5)) = qm(Q) (where q = | det In 
particular, 

qm{T) = m(A(T)) = m I [j (T + d) j (by property 2 of T) 

\d£D ) 

< m ( T + d ) = Yl m ( T ) = qm{T). 

deD deD 

Additionally, property 2 implies that A k+1 (T) = \J deD (A k (T) + A k d) for all k > 0. 
Then 

g fe+1 m(T) < m(A k (T) + A k d) = ^ q k m(T) = q k+1 m{T). 

deD deD 

Thus 

m ((^ fc (T) + A fc dj) f|(^ fe ( T ) + A%)) = 

for all distinct di, dj G D. 

Since 6 D, A k+1 (T) D A k {T) for all k > 0. So 

(4*+i(T) + A fc+1 d) D (A fe (T) + A k+1 d) 

for all d £ D. Now consider T + x for any i £ Z". By hypothesis, 

N-l 

x = A N d N + A ° d 3 

3=0 

for some TV > 0, with d N ^ G D and the G £>. Then 

JV-l -1 

T + x = {y G R" : y = A N d N + ^ A j dj + ^ Aid J> a11 d : e 13 } 

3=0 3=-oo 

C J 4 JV " _1 (T) + A^dw c A W (T) + A^dAT, 

and 

(Tf|(T + a:)) C (A W (T) f](A N (T) + A N d N )). 
Then d^ / OeD implies that 

= m (a n (T) f]{A N {T) + A N d N )^j > m (rf)(T + *)) . 

□ 
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Figure 1 . The tile T for the twin dragon matrix. 



If the set T = T(A, D) tiles K™, it is not necessarily true that A yields a radix 
representation. For example, in the case where A = 2, we can find a radix rep- 
resentation for all nonnegative integers with the digit set D — {0, 1}, o r for all 
non-positive integers with the digit set D = {0,-1}, but we cannot represent all 
integers with a radix representation using any digit set [TO]. Yet T = [0, 1] (with 
digit set D = {0, 1}), does tile R under translation by Z. Similarly, if A is the twin 
dragon matrix [12] . 

11 

-1 1 



A = 



then A is a dilation matrix. A digit set for A is 





D = {d = 







di = 



}■ 



and the set T generated by the twin dragon matrix with this digit set tiles R 2 
under translation by Z n [12] , yet A does not yield a radix representation of Z 71 (for 
example, the vector 


-1 

does not have a radix representation) [TJ. 

In both of the examples, the origin is on the boundary of the tile T, so that 
U£L A fe (T) C R™. We claim in Theorem [0J that this must be the case for all 
dilation matrices A which give tiles T but do not yield radix representations. We 
first introduce a technical lemma. 

Lemma 8. Let A be a dilation matrix, and let D be a digit set for A. If T° 
then there exists an increasing subsequence {£j}j>i of the positive integers and a 
sequence of vectors {Q : Q € A~ lj (Z n )}j>i converging to such that Q $ T for 
all j > 1. 

Proof. In the proof of Lemma [7] we noted that the set T is compact. Thus for any 
uj £ R™ such that u> £ T, d(u>,T) > 0, and there exists an open ball centered at u>, 
B u , such that T = 0. 

Let {yj}j>i be a sequence of vectors in R™ converging to with yj £ T for all 



j > 1. Set ej = \\yj\\p, and notice that €j > d(yj,T) > 0, with lhxij 



0. 
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Set n = %^ and 

. (rj-i d(y j ,T)\ 
ri =nun| — 1 

for j > 2. Then {r^}j->i is a decreasing sequence of positive numbers, lim : j_ > . 00 rj — 
0, and the open balls Bj of radius rj centered at the vectors yj satisfy BjC\T = fy 
for all j geql. By construction, if y* is any point in the ball Bj for each j > 1, then 
linij-^oc y* = and y*j^T for all j > 1. 

In particular, there exists an increasing subsequence {(j}j>i of the positive in- 
tegers such that A~ lj is a fine enough lattice to ensure that A~ £j (Z™) D Bj ^ 
for each j > 1. We may choose some (*,- £ /l~ £j (Z n ) n for each j > 1. By 
construction, (j T for each j > 1, and the sequence {(j}j>i converges to 0. □ 

Lemma 9. Let A be a dilation matrix, and let D be a digit set for A. If A yields a 
radix representation with digit set D, then is in the interior of the set T generated 
by A and D . 

Proof. We prove the desired result by contradiction. Assume that ^ T°. By 
Lemma [51 let {Q}j>i be a sequence of vectors in R n converging to such that 
C,j\nA~ l] {JL n ) (for some increasing subsequence {£j}j>i of the positive integers) 
but £j £" T for all j > 1. We may write Cj = A~ lj Xj for some Xj £ Z™ for each 
j. Since A yields a radix representation, there exists an integer Nj for each Xj and 

digits . . . , d$. £ D such that 

Xj = Y,A i d<?\ 

i=0 

Thus 

AT, 
i=0 

with % £ Z" and ^ £ T. 

Since T is compact, ||£jl|; 2 is bounded above by b for some b > 0. Then, since 
converges to 0, is also bounded above for sufficiently large j. We use 

the rough estimate that there exists an intcg er M > 1 such that for all j > M, 
\\kj\\[2 < 26. Thus for j > M, the integer vectors kj all belong to a finite subset 
of Z™. This implies that there exists an integer N > such that Nj < N for all 
j > M. For j > max {M, TV}, N - £j < N - j < 0, and thus ( 3 £ T for all 
sufficiently large j. This contradicts the choice of Q, thus our assumption that 
£■ T° must be false. □ 

We have shown that if a dilation matrix A yields a radix representation with digit 
set D, then the set T = T(A, D) tiles l n under translation by Z", and £ T°. 
Next we prove the converse. 

Lemma 10. Let A be a dilation matrix, and let D be a digit set for A. If the set 
T = T(A,D) tiles R™ under translation by Z™ and ifO £ T° , then A yields a radix 
representat ion with digit set D. 

Proof. Following the notation of [6], set 

k-l 

Dam ■= {x £ Z" : x = ^ A ] dj,d £ D}, 

3=0 
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the set of vectors that can be expressed with a radix representation of length less 
than or equal to k. Note that Da.i = D. By construction, A(T) = UdeD{T + d). 
Thus 

A k {T)= |J (T + x). 

xeD A .k 

By the tiling hypothesis, (T° + x) n (T° + y) = for all distinct x,y e1 n . Thus 
(T° + y) n A fc (r°) = for all y e Z" with y £ £) Afe , and 

D A , k D(A k (T°)f)Z n ). 

Let B be an open ball centered at the origin such that B C T°. Then 

D A ^ k D(A k (B)nZ n ). 

The sets A k+1 (B) are expanding, with U fc >o^ fc+1 ( B ) = R "- Tll us 

U^AfeS U(^( s ) nz ")= z "- 

fc>0 fc>0 

The opposite containment is true as well, since Da,u C Z™ for each fc. Thus 
U fc > -DA,fe = Z n - □ 

We have now completed the proof of Theorem [BJ 

Recall that \i is the smallest singular value of the dilation matrix A, and that F 
is our canonical fundamental domain of Z n , F = [— |, -|) n . Combining Theorems [3] 
and [71 we also have t he following corollary. 

Corollary 11. Let A be a dilation matrix, and let D = A(F) n Z™. If fi > 2 then 
the set 

-l 

T={ieM™ : x = ^ A j dj, dj € £>} 

j=-oo 

fifes M" under translation by Z" . 

3. Haar-Like Wavelets 

Sclf-affine tiles allow us to construct multivariable wavelet sets associated with 
multiresolution analyses. We review some basic definitions from wavelet theory 
here. 

Definition 12. A multiresolution analysis (MRA) associated with a dilation ma- 
trix A is a nested sequences of subspaces • • • C VLi C Vo C V\ C • • • of L 2 (R n ) 
satisfying: [12] 

(1) U^ = L 2 (M»); 

(2) nj&Vj = {0}; 

(3) f(x) G Vj if and only if f(Ax) e v j+i for all j e Z; 

(4) f{x) £ Vo if and only if f(x - k) E V for all k e Z"; and 

(5) there exists a function </>(x) 6 Vo, called a scaling function, such that 

{^(a; - k) : fee Z"} 
is a complete orthonormal basis for Vq. 
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The existence of multiresolution analyses in dimension n > 1 has been stud- 
ied by a number of authors. In [3], Grochenig and Madych showed that if cf> — 
(rn(Q)) 1 ^ 2 XQ is a scaling function for a multiresolution analysis, then Q mus t be 
an affine image of a self-affine tiling of K" under translation by Z". They showed 
also that if T is a set of the form 

-l 

r = {i£R": x = ^ ,l"',/,.</, 6 £>} 

j=—oo 

with A a dilation matrix and D a digit set for A, then = \t is the scaling function 
for a multiresolution analysis (note that \T\ = 1). Lagarias and Wang noted that all 
self-affine tiles T which tile R" under translatio n by Z™ must be of this form [8]. A 
scaling function that is the characteristic function of some measurable set is called 
a Haar-like scaling function (after the Haar scaling function, which is X[o,i])- 

Lagarias and Wang studied necessary conditions for sets of the form T(A, D) to 
tile W 1 under translation by Z™ in a series of papers ([5], [TJ, [5], [5]). Much of their 
work in these papers concerned t he question of when a set T(A, D) tiles M. n under 
translation by a sublattice of Z n . In [6] and [TJ, they showed that if D is a complete 
set of coset representatives of Z n /A(Z n ), then T(A, D) is a self-affine tile of W 1 
under translation by some sublattice of Z™ . They also studied some properties of 
the tiling set T{A,D). He and Lau g] studied sets T(A, D) which tile R n under 
translation by a sublattice of Z™ as well; in particular, they looked at possible digit 
sets D. 

In order for \t to be a scaling function for a multiresolution analysis, however, 
we need T to tile M™ under translation by all of Z n . Lagarias and Wang gave some 
necessary conditions for a dilation matrix A to yield a Haar-like scaling function 
in [5] and [5]. They showed that all dilation matrices in dimensions n — 2 and 3 
yield Haar-like scaling functions. Our results below give a sufficient condition for a 
dilation matrix A to yield a Haar-like scaling function, in any dimension. Note that 
for dilation matrices that yield a Haar-like scaling function, Strichartz has shown 
that multiresolution analyses and associated wavelet bases with arbitrary regularity 
can be constructed [TT] . 

The results of the previous section imply the following two theorems. 

Theorem 13. Let A be a dilation matrix, and let D be a digit set for A such that A 
yields a radix representation for Z™ with digit set D. Let T be the set depending on 
A and D defined above. Then (f> — \T is the scaling function for a multiresolution 
analysis. In particular, if A satisfies [i > 2 and if D is the set D — A(F) flZ" with 
F = [— i, i) n ; then cj) = \t is the scaling function for a multiresolution analysis. 

Theorem 14. Let A be a dilation matrix. Then there exists a positive integer f3 > 1 
such that for all integers k > fj there exists a multiresolution analysis associated 
with the dilation matrix A k . 

Thus Haar-like scaling functions and associated MRAs exist for a large class of 
dilation matrices. 
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